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Abstract 

Motivated by the structure of one-loop vacuum polarization effects 
in curved spacetime we discuss a non-minimal extension of the Einstein- 
Maxwell equations. This formalism is applied to Bianchi I models with 
magnetic field. We obtain several exact solutions of the non-minimal 
system including those which describe an isotropization process. We show 
that there are inflationary solutions in which the cosmological constant 
is determined by the non-minimal coupling parameters. Furthermore, we 
find an isotropic de Sitter solution characterized by a "screening" of the 
magnetic field as a consequence of the non-minimal coupling. 



PACS numbers: 04.40.Nr, 98.80.Jk 



1 Introduction 

The Einstein-Maxwell theory has been a subject of investigations since long pQ. 
As far as cosmological implications are concerned, the possible role of a primor- 
dial magnetic field has attracted particular interest. Studying the impact of 
a magnetic field on the dynamical evaluation requires anisotropic cosmological 
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models. A general discussion of this type of models with magnetic field and 
references to early activities along this line may be found in By now there 
are strong limits on the current magnitude of such a field 3 which seems to 
render magnetic fields on cosmological scales unimportant at the present stage 
of the cosmic evolution. This does not imply, however, that a magnetic field did 
not influence the dynamics of the early Universe. Moreover, quantum effects are 
expected to become relevant at early cosmic stages. Quantum electrodynamical 
consideration show that vacuum polarization effects in curved spacetime give 
rise to non- minimal modifications of the (minimal) Einstein-Maxwell Lagrangian 
0]. The investigation of a non-minimal coupling of gravity with electromag- 
netic fields was initiated by Prasanna 5 . Prasanna introduced the additional 
invariant R lkmn Fif.F mn (R lkmn is the Riemann tensor, Fik is the Maxwell ten- 
sor) into the Lagrangian for the gravito-electromagnetic system and obtained 
a non-minimal one-parameter modification of the Einstein-Maxwell equations 
|S] . Novello and Salim p] included the (gauge-dependent) terms RA k Ak and 
R lk AiAk in the Lagrangian {Ak is the electromagnetic potential four-vector, 
R lk is the Ricci tensor and R is the curvature scalar) . A qualitatively new step 
has then been made by Drummond and Hathrell 0] by calculating quantum- 
clectrodynamical one-loop corrections in curved spacetime. The Lagrangian of 
such a theory contains the three U(l) gauge-invariant scalars R lkmn F i kF mn , 
R lk g mn F im Fk n and RF mn F mn with coefficients proportional to the square of 
the Compton wavelength of the electron. Subsequently, a non-minimal coupling 
of gravity and electromagnetism has been discussed by a number of authors in 
different settings 0131101 El El El HI EH Gl DUGS!. Non-minimally extended 
theories were used as a framework to discuss potential limitations of the equiva- 
lence principle |19l 1201 I2T] ). A further quantum electrodynamical motivation of 
the use of the generalized Maxwell equations can be found in |23 EH1 12] • The 
effect of birefringence induced by curvature, first discussed in jj] , and some of its 
consequences for electrodynamic systems have been investigated for the pp-wave 
background in j^SJ [2H1 H3 1281 . A curvature force has been introduced to 
describe the accelerated expansion of the universe |3l)l l-'Slj . Non-minimal inter- 
actions in which torsion is coupled to the electromagnetic field were studied in 
123123 (see also for a review). Finally, we mention a mediated non-minimal 
coupling in which the scalar Higgs field <f> is coupled to gravity via a £,(j) 2 R term 
and to a Yang-Mills potential Ak by (j) 2 AkA k |3"5] , 

Most of the investigations so far were devoted to the analysis of the non- 
minimally modified Maxwell equations on a given background. With the excep- 
tions |5] and |15) the impact of the non-minimal coupling on the gravitational 
field has been outside the focus of interest. However, one may expect that 
the rich structure of such type of theories gives rise to novel features in the 
gravitational dynamics as well. Our purpose here is to demonstrate this aspect 
for a specific class of non-minimal interactions in Bianchi I cosmological mod- 
els. These interactions are modelled according to the already mentioned general 
structure obtained in as the result of quantum electrodynamical calculations. 
While the coupling constants are fixed in ^ , they will be considered as arbitrary, 
constant parameters in our analysis (cf. ^SJ)- On this basis we shall first obtain 
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the general set of equations for a non-minimally extended Einstein-Maxwell sys- 
tem with linear electrodynamics (For a non-linear, non-minimal extension of the 
Einstein-Maxwell theory see [SHI)- Then we specify this set to the homogenous 
but anisotropic case of Bianchi I cosmological models with magnetic field and 
a matter component with generally anisotropic equations of state. We obtain 
simple exact solutions for several choices of the non-minimal coupling param- 
eters. Among them are solutions with axial symmetry which isotropize in the 
long time limit. As a specific feature of inflationary solutions of the non-minimal 
theory we find direct relations between a cosmological constant and the non- 
minimal coupling parameters. There exists an isotropic de Sitter solution with 
this property as well for which the corresponding set of parameters makes the 
gravitational dynamics independent of the (non-vanishing) magnetic field. 

The paper is organized as follows. In section[2]we present the general formal- 
ism of the non-minimally extended set of gravito-electromagnetic field equations, 
based on 0] . The (linear) electrodynamical field equations are obtained an dis- 
cussed in section [3J Section 0] is devoted to the gravitational field equations in 
the general case. In section [3] the material of the previous sections is applied to 
the Bianchi I geometry. Several particular models are then studied in section 
Section [7] provides a summary of the paper. 

2 Non-minimal coupling of gravity and electro- 
magnetism 

2.1 General formalism 

A non-minimal extension of the Einstein-Maxwell theory can be derived from 
the action functional 



S\s,A 




with the Lagrangian 

D _j_ O A 1 1 

t ~*~ if I _ jp jpmn _i_ _ijifcmn jp jp (n\ 

^ ~r ^matter \ r . r mn r T" ~ /V r ik r mn • \^ ) 

Here, A is a cosmological constant, g is the determinant of the metric tensor 
gik, the constant k is equal to k = ^f- where G is Newtons gravitational con- 
stant. The quantity £ ma ttcr is the Lagrangian of neutral matter and is the 
Maxwell tensor = Vi^4fe — \7 k A i: where Vj denotes the covariant derivative 
and Ai is the four potential. The last term describes a (U(l) gauge invariant) 
non-minimal coupling between gravity and electromagnetism, mediated through 
the tensor 

n ikmn = ^ g im g kn_ g in g km^ R 

+^(R im g kn -R in g km +R kn g im -R km g™)+q 3 R ikmn , (3) 
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where q%, q 2 and q$ are phenomenological coupling constants with the dimen- 
sion [length] 2 . This structure of the coupling is motivated by quantum electro- 
dynamical calculations of vacuum polarization effects in curved spacetime by 
Drummond and Hathrell While qy, q 2 and q^ have definite values in 0], we 
assume them to be arbitrary constant parameters in our analysis. The choice 
© is a generalization of previous non- minimal modifications of Maxwell's 
theory. The case q\ = q 2 = was investigated in |SJE]. For q 2 = qz = one 
obtains a model considered in |15) . 

The tensor TZ lkmn has the same symmetry properties as the Riemann tensor 
R lkmn , Contraction yields 

g kn Tl ikmn = R im ( q2 + q,) + ^Rg" n (3q 1+ q 2 ), 

gkn9imK ikmn = R(6qi+3q 2 + q 3 ). (4) 



The case of a vanishing trace which will be of interest in later applications is 
characterized by 

g kn g tm n tkmn ^ => 6 qi +3q 2 + q 3 =0. (5) 



3 Electrodynamic equations 

The equations of non-minimal electrodynamics are obtained by varying the ac- 
tion functional with the Lagrangian J2J with respect to the four-potential Aj of 
the electromagnetic field. They are of the standard form 

V k H ik = , V k F* lk = , (6) 

where the induction tensor H lk is given by 

H ik _ F ik + Tz ikmn F mn , (7) 

This relation has the structure of a constitutive law in which fZ' lkTnn plays the 
role of a susceptibility tensor. 

3.1 Non-minimal constitutive equations 

The linear constitutive equation Q has the standard form H lk = C %kmn F mn 
123 EH1 EH with a "material" tensor 

^jikmn ^_^g'^ rn g^ n gitigkm^ _|_ j^ikmn 
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This tensor describes the linear electromagnetic response of the system, which 
may also be characterized by the dielectric and magnetic permeabilities, as well 
as by possible magneto-electric effects [37| EHI HUD- C lkmn can uniquely be 
decomposed with respect to the four velocity U l (normalized by U l Ui = 1) of 
the medium: 

^jikmn jj_ ^imjj-kjj-n ^injjkjjin _|_ ^knjjijjm ^kmjjijjnl^ 

-\rf H ^- X )ur, mns 

_ 1 \jfkl(jjm Vi n _ un v m-j + ^mn^ji^k _ U^*.)] . (9) 

Here e lm and {n~ 1 ) pq are the dielectric and magnetic permeability tensors, re- 
spectively, and is p m is a tensor of magneto-electric coefficients. These quantities 
arc defined as 

e im = 2C lkmn U k U n , {^% q = -^kC tkmn Vmnq , (10) 

v™ = T)pikC ikmn U n = U k C mkln mnp . (11) 

The dot denotes the position of the second index when lowered. The tensors 
r) mn i and rf kl are anti-symmetric and orthogonal to U l , 

Vmnl = trnnlsU 8 , rf U EE e Ms U s . (12) 

They satisfy the identity 

- V lkp V mnp = C ns UiU s = A) n A k n - A^At , (13) 

where 8J^ na is the generalized 6- indices 8— Kronecker tensor. The spatial pro- 
jection tensor A lk is defined by 

A ik = g ik _ jjijjk _ ^ 

Contracting equation l(T3|l yields 

\v m mi m = -€ s UiU s = -A* m . (15) 

The tensors and (/i -1 ),:^ are symmetric, but v\ k is generally non-symmetric. 
All these tensors are orthogonal to U l , 

e lk U k = 0, {iT 1 ) ik U k = 0, v k .U l =0 = v k .U k . (16) 

Use of (|Q in (fTTTfl and |(TJ| provides us with 

£ vm = ^im + 2K tkmn U k U n , (17) 
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v v m = vpi k n ikmn u n . (19) 

The non-minimal coupling of gravitational and electromagnetic fields may be 
interpreted as a change of the dielectric and magnetic properties of the vacuum, 
including a specific magnetoelectric interaction. The vacuum acquires prop- 
erties of a quasi-medium under the influence of a non- vanishing tensor TZ lkmn . 
The analogy between non- minimally extended electrodynamics and macroscopic 
media was pointed out, e.g., in [22]. This analogy may be completed by intro- 
ducing the electric induction D\ the magnetic field H l , the electric field E l and 
the magnetic induction B l , [37]: 

D* = e im E m -B l v l \, H i = u i m E m + { l x- l ) im B m . (20) 

The vectors D\ H\ E i and B i are defined by gU]: 

D i = H ik Uk ^ Rl = H *ik Uk ( E i = ^fc^ ^ fli = F *ik Uk _ (- 21 ) 

They are orthogonal to the velocity four- vector U l : 

D l U t = = E l U, , H l U t = = S l C/ 4 , (22) 
and form the basis for the decomposition of the tensors F mn and H mn : 

Emn — E m U n E^-JJfn Tj mn lB , H mn — D m ,U n D n U m Tj mn lH . (23) 



4 Gravitational field equations 

The gravitational field equations are obtained by varying the action Q with 
the Lagrangian ^ with respect to the metric tensor. They can be written in 
the standard form 

R ik ~lRg tk = Ag lk + kT^ S) , (24) 

where 

4f 3 = Tt 1 ^ + 4° } + rf } + ^ + fl»T#> ■ (25) 
The stress-energy tensor of the matter 2^™ attcr ) mav De decomposed according 



^(matter) = + ^g) ^ + ^ ^ + ^ ^ (2g) 



to 

where is the matter energy density scalar, is the symmetric (generally 
anisotropic) pressure tensor, orthogonal to the velocity four- vector {P ik U k = 0), 
and l[ q ^ is the energy-flux four- vector, orthogonal to the four velocity (I^ q 'U i — 
0). By T$ we denote the usual stress-energy tensor of the electromagnetic 
field, 

T ik = -^9ikF mn F mn - F m F k n . (27) 
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The contributions from the non-minimal interaction are 



T 



(i) _ 



R T- 



(0) 



1 



Hk" mn-L 



(28) 



T, 



-i?"^ im F fc „--V'V ; (^ m F fe ") 

+ iv ; [V 4 (F fc „F / ")+V fc ( J F mJ F lta )] , 



(29) 



and 



1 *3 

T>(3) ~ r> mills rp TP T?ls ( T?' n T? i J?' n £? \ 

1 ik — ~^9ikM * i mn- t ls — -7-F i ^knls T ^ fc Minis) 

-iv m V n (i^ fe ™+F fe "i^). 



(30) 



The tensor Tj^ is proportional to the corresponding term in ^5], the part 
-Vfc reproduces the stress-energy tensor of The tensor T^ fc ; is new. All 

three terms are supposed to contribute to the total stress-energy tensor in the 
following. In contrast to the traceless electromagnetic stress-energy tensor T t 



the tensors T-£' , T^ 1 and TrP have non- vanishing traces: 



(0) 

ik 



y i ik 



= -9i 



ikrp{2) 

9 1 ik - ~12 



ikrp(3) 



"?3 



^RF mn F mn + ^V k V k (F mn F mn ) 
R mn F k m F kn + iv fe V fe (F mn F™ n ) 



nmnls tp tj i Y7"iV7 / Tjkn t? \ 



(31) 
(32) 
(33) 



Non-vanishing traces of effective stress energy tensors are also features of non- 
linear electrodynamic models (see, e.g., |41|). 

The effective stress-energy tensor (|2*5|) in eq (|2"^ has to be divergence-free, 

i.e. 

krp(eS) _ 



ik 



. (34) 

We assume the stress-energy tensor of the matter 7 l .(™ attcr ) to be conserved 
separately, i.e., v fe y^™ attor ) = o. The remaining part of the effective stress- 
energy tensor is then automatically conserved if F ik is a solution of Maxwell's 
equations. In order to check this fact directly, one has to use the Maxwell 
equations 10 with (JJJ, the Bianchi identities, the symmetry properties of the 
Riemann tensor and the commutation rules for the covariant derivatives. This 
procedure is analogous to the one described in |lfi| and we omit it. 
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5 Anisotropic cosmological models 

5.1 Metric structure 

We consider now the Bianchi I cosmological model with the line element ^ |5] 

ds 2 = dt 2 - a 2 (t) (dx 1 ) 2 - b 2 (t) (dx 2 ) 2 - c 2 (t) {dx 3 ) 2 . (35) 

Due to the symmetry of the metric only six components of the Riemann tensor 
arc different from zero: 



nOl _ _ O p02 _ _ ^ r>03 _ ^ 
-"-■■01 ~~ ' rL -02 — T 7 rt -03 — ~ 



pU « 5 P 13 «C t-,23 fo C /o fi N 

K -12 — 7 J H -13- ■ W -23- _ 77- W 

a o a c c 

5.2 Exact solution of Maxwell's equations 

Since the susceptibility tensor TZ lkmn has no non-vanishing components with 
only one index zero (TZ a ^ = 0, Greek indices denote spatial coordinates), 
it follows from relation (|19f) that all magneto-electric coefficients vanish. The 
gravitational field (|35[1 does not mix pure electric and pure magnetic fields. 
From relations l|17|l and l |18|l we find the dielectric and magnetic permeability 
tensors 

Let us consider now a magnetic field directed along the Oz axis (which is also 
the direction of the shear eigenvector 02). The symmetry of the problem then 
fixes Fik and H lk to be of the following structure: 

F ik = (SIS 2 - S 2 5\)F 12 , (38) 

H lk = (g a g k2 - g i2 g hl )[l + QiR + <fc(ie} + R%) + <feie? 12 ]fi 2 • (39) 

The second set of the equations © yields F 12 = const. The first set of eqs. © is 
identically satisfied since all the components H t0 are equal to zero. Introducing 
the scalar value of the magnetic field B(t) by 

B 2 (t) = ^F lk F lk = F 12 F 12 , (40) 



we reproduce the result [cf. [2])] 



i.e., B(t)a(t)b(t) = const. 
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5.3 Einstein's equations 

In a co- moving frame with U l = 6$, the field equations l|24|l reduce to the 
following system: 



a b at be 

-T + -- + T- 

a b a c be 



A + nW+ -nB 2 {t) 



+ K B 2 (t) { qi 




ab c I a b 



92 



a b 
a b 



a b 
a b 



(42) 



c be 
c be 



A-kP, 



(i) 



-nB 2 (t) 



-<T2 



, a c 



4-- +4 

ab a c 




b c 

ab c I a b\ 
-2-T+2- -+- 
ab c \ a b j 



53 



Hi 



be 
be 



(43) 



c a c 
c a c 



A-kP, 



(2) 



nB 2 (t) 



+ K B 2 {t) { Ql 



b a 
4- + 3- 

o a 



+92 



a 6 




,06 ,6c „ a c 
4_ +4 +3__ 

a o be a c 



ab eld b 
2-T+2- -+t 
a o c \ a o 



a__ 2 la 



a c 
a c 



, (44) 
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r/ b 
a b 



A — kPi 



+ K B 2 {t) { gi 




(45) 



Here, P(i), P12) and P(3) are the eigenvalues of the anisotropic pressure tensor 
Pik- Summing up eqs. I|42|l - H45|) we obtain the trace equation: 



ft 



o. b 



a c c 



a b c a b a c be 

->2 



4A + K(W-P il} -P {2} -P {3) ) 



-kB (t) { (8 0l 



1i + 03 



-2(6gi + 3q 2 + ffi) 




2(5?i + 2g 2 )-^ 
a 



(8gi+4g 2 +g 3 )- ( - + - 
c \ a b 



(46) 



Differentiating (|42|l and using (|43H - 145|l leads to the conservation law for the 
matter: 



W 



a b c\ a 
- + - + -}W+ -P (1) 
a c a 



-P, 



(2) 



-P, 



(3) 



0. 



(47) 



It follows that + qi T^> +q 2 T^> +q 3 T^> is conserved as well. The equations 
l|42l) - (|45|l represent a modified (compared to minimal coupling) dynamical 
system, since the non-minimal terms contribute to the coefficients before the 
second order derivatives. Its complete analysis should be the subject of future 
investigations along the lines described in |43| . 



^(1) 



(3) 



6 Particular models 

6.1 Quasi-one-dimensional solutions with pure magnetic 
field 

As the first application we consider a model without matter (W = 0, Pit = 0) 
and with constant values of two of the metric functions: a(t) = oq and b(t) = 6q. 
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Only c(t) is assumed to vary. Then B 2 (t) = const (see, Q4H ) and the equations 
JUJ-gSJ reduce to 



= A 



-nB 2 {t ). 



- c [l-q lK B 2 (t )] = A-^B 2 (t ) 



(48) 



Only one of the non-minimal parameters enters the dynamics. According to 
the first equation this model requires a negative cosmological constant which 
exactly compensates the magnetic field term, i.e., kB 2 = — 2A. The equation 
for c(i) is 

KB 2 (t n ) 



%t) + c(t) 



1 - qi KB 2 (t ) 



0. 



(49) 



For B 2 = A = we obtain the vacuum solution c oc t which is the degenerate 
case of a Kasner solution, equivalent to the flat spacetime Milne universe (cf. 
[2]). For non-vanishing kB 2 — — 2A we have three cases. 

First case: q\nB 2 (to) < 1. 

This condition includes the case q\ = 0. The solution of eq. I|49|l is oscillatory: 



where 



c(t) — c(to) cos v(t — t ) + £-^L s i n v (i _ £ Q ) 



k£? 2 (£ ) 



l- giK B 2 (t ) 
The zeros of c(t) denote singularities of the model 



(50) 
(51) 



Second case: qiK,B 2 (to) = 1. 

The equations l|48l) contradict each other. This case is incompatible with B(to) ^ 
0. 

Third case: q\KB 2 (to) > 1. 

This condition requires qi to be positive. The solution of (|4^|l is 



c(t) = c(t ) cosh — t ) + £_^2_ s i nn ^(t — 1 ) , 



(52) 



where 



2 _ K g 2 fa) /„x 

" U>^ 2 (*o)-lJ ' ( j 

For t >> to the function c(t) behaves as c(t) oc e Mt . The model is non-singular 
when | C y^ ) 1 > 1 - If | °Mtl) 1 < ■*-! tnen there exists a singularity at t* , given by 
tanh^(r - t ) = and c(f*) = 0. 



c(to) 
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6.2 Quasi-two-dimensional solutions with magnetic field 
and matter 



A second model with B(t) =const is obtained for a(t)b(t) =const. We may 
write a(t) — a(t )E(t) and b(t) = 6(to)i? _1 (i) with E(t ) — 1 where to is some 
reference time. If we additionally assume c(t) — const, the dynamics is restricted 
to the x 1 Ox 2 plane. The equations (|42[) -H45 [1 reduce to the system: 

-l(^) = A + nW+^B 2 , -i| + 2 f| N ) = A- K P {1) -±kB 2 , (54) 



= A-kP (3) + ±kB 2 , L~ + 2(l-L)(~") =K-kP [2) - 1 -kB\ 



(55) 



where 

L=l + (qi- q 3 )nB 2 = const . (56) 
For L^O the system is equivalent to 

2L (f\ =-«(W + P (3 )) > 2L (§) = <P(i) - P(2)) , (57) 



P(3) —W + B 2 -\ , P (1) + P {2) = 2W H . (58) 

The minimally coupled case q± = q 2 = q 3 = corresponds to L = 1. It is 
dynamically indistinguishable from a non-minimal configuration with q\ = q 3 
and arbitrary q 2 . The case L = 1 requires W + Pr^ < 0. This can only be 
achieved by a sufficiently negative cosmological constant, 2W + B 2 + ^ < 
which also implies that Pm + P( 2 ) < 0. 

6.2.1 Case L = 

For q 3 — q\ = — the quantity E(t) is constant, i.e., the universe is static. The 
matter distribution is characterized by constant quantities as well: 

P {x) = P{2) = ±- l -B\ W = -P m = ~-\B\ (59) 

The energy density W is positive if A < —^kB 2 . All the pressure eigenvalues 
become negative in this case. 
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6.2.2 Case L^O, P {1) = P { 



(2) 



Here we obtain 



E(t) = e Ho{t - to) , 



K 



P(3) — ~LHq H h -B 2 



P(i) — P(2) — -LH + — - -B 



A 1 



(60) 
(61) 



where i?o is an arbitrary integration constant. For any L > a sufficiently 
negative cosmological constant is required for the energy density W to be posi- 
tive. After re-parametrization of the coordinates and the time the metric takes 
a form 

ds 2 = dt 2 - (e^dx 2 + e-^dy 2 ) - dz 2 . (62) 



6.2.3 Ultrarelativistic matter with L^O, P(i) = P( 2 ) 

If the matter is ultrarelativistic, i.e., TK mattc ^ = and, consequently, W 
P(i) + P<2) + P<3), onc obtains 



„2 2A 



W 



-Ilh 2 -\b 2 

2 2 



(63) 



P, 



(3) 



-LH 2 



-B 



P(D = P i 



1 



(2) 



LH 2 



-B 



(64) 



The energy density W is positive for LHq < —^B 2 , which requires the constant 
L to be negative, i.e., q% — qi > In this case the longitudinal pressure 

P( 3 ) is positive as well. The transversal pressure P<\\ = P/ 2 \ is positive for 



I 2 < -B 2 . 
constant again. 



LHq < — B z . The condition L < necessarily implies a negative cosmological 



6.3 Axial symmetry: a(t) = b(t), = P( 2 ) = -P( tr ) 

All the metric functions are assumed to be time dependent now. Einstein's 
equations reduce to the following system of three equations for two unknown 
functions: 



- j 2 + 2— - = A + kW + kB 2 (£) i tt + (7gi + 4g 2 + q 3 ) ( * V + 2q^- 
a I ac 2 \a ) ac 



(65) 



- + - + --=A-KP (tr) + K B 2 (t) { --+?!-+ 7gi+4 g2 +g 3 
a c a c v ' | 2 c 



a „ / a \ a c 
—2 - +-- 

a \ a / a c 



2 ~a + Ui) = A -«- P (3)+«- b2 W | l + 2 ^l- (l) (13gi +4g 2 + g 3 ) I - (67) 
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The matter conservation law takes the form 

W + 2-(W + P (tl) ) + -(W + P (3) ) =0. (68) 



Formally, the function c(t) does not appear in equation l|67|) . However, in general 
there is a coupling to c(t) via the pressure P( 3 ) through the conservation law 
<|6l31) . There are two simple cases for which eq. ljrJ7|) decouples from eqs. l|63|l 
and H66f) . The first one is P( 3 ) = 0, i.e., a vanishing longitudinal pressure, the 
second one is W+P^ = 0, i.e., a vacuum type behavior in longitudinal direction 
together with a transversal equation of state P( tr ) = P( tr ) (W) . In the following 
we consider both cases separately. 

6.3.1 "Longitudinal dust": P (3) =0 

With the substitution 



a(t) = a y/Z(a) , cr = , B (t) = — — , (69) 

Zqi a 4 

where M 2 =const and qi ^ 0, equation (|67|) is transformed into the first order 
equation 

n qi M 2 \ dZ(a) , ,Z / kM 2 \ , „„ 

^VJx^^^rr^J a ■ (70) 

The solution of eq. I|7U|I can be represented as quadrature: 

2<T + 1 

^(«) = H a(ao)a 



2 



(71) 



a 4 - nq\M 2 

a 4 - KgiAf 2 )"^ / dxa; 1 " 2 ' 7 (x 4 - Kqi M 2 )^ [Ax 4 + \nM 2 
Jan V 2 



where H a (a Q ) = a^V^M with ao = a(i ) and H a (t) = f = £T t (i) = £ 
is the expansion rate in the x 1 and a; 2 directions. Due to the expressions 

(a; 4 - nqiM 2 )^' and (a 4 - KqiM 2 )^~ in Eq. (JHJ, the result of the inte- 
gration is sensitive to the sign of q\. For qi > the term (cc 4 — /cgiM 2 ) has 
two real zeros cci^ = ±(KgiM 2 )3. If at least one zero belongs to the interval 
(ao,a(t)), the integral in eq. H3} diverges for 2 ?f^ < -1, i.e., for 2a + 1 < 0. 
If q% < 0, such a singularity does not appear. The asymptotic behavior of the 
function Z(a) for a — > oo is 

Z(o -» oo) = ^a 2 ~ 2CT . (72) 
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Thus, asymptotically all the models yield 

- = J- = H a = const , a(t) = a^e"^-^ , (73) 
a V o 

i.e., a de Sitter type expansion for a(t), independent of the parameter a. More- 
over, one can check directly that the solution l|73|) is an exact solution of the 
equation (|67|1 with 

H- 2 = 2(11^+4(72 + 93), (74) 

where H91 + 4g2 + 93 + 1 0. It is remarkable, that the constant expansion rate 
H a is determined both by A via eq. I|73|) and by the non-minimal coupling 
parameters via eq. Ij74(l . The relation Q74JI has no counterpart in the minimal 
theory. Combining (|73|) and l|74|) allows us to establish the following relation 
between A and the coupling parameters: 

3 

A = 2(llgi+4g 2 + g3) ' (?5) 

The cosmological constant is expressed in terms of quantities which are sup- 
posed to be the result of quantum field theoretical calculations. In quantum 

electrodynamics the parameters gi, 92, and 93 are A, 91 = — 92 = — 13gi, 
93 = 291, where a is the fine structure constant, i.e., they are proportional to 
the square of the Compton wavelength A e of the electron. This would give rise 
to a value Aqed = fg^i 7 °^ tne cosm ological constant. While one does not 
expect a quantum electrodynamical length to set the scale for an early de Sitter 
stage (recall that we assume 91, 92, and 93 to be free parameters) this result 
may nevertheless indicate a potential relevance of non-minimal interactions for 
an early inflationary dynamics. 

With a transversal equation of state P( tr )(i) = (7 — l)W(t) and with l|73|l 
the conservation law (|68f) yields 

W {t) = Witof^-e- 211 ^-^ . (76) 
The equation Q65[l for c(t) can be rewritten in the form 



where 



ffz {[1 - az]-<Y{z)} = -[I a Z ]-^^>z^ , ( 77 ) 



-\a(t)J ' - 4 ' 891^ ' -a(t)c(t )- [<S) 

For a 7^ 1 the solution for c(t) with the initial value c(to) has the following 
explicit form 



c(t) = c(t )e^( t -*«) 



1 _ ae -4ff.(t-to) 

4H a (t-t n 



(i-«r e -^/ ■ — r ,1+l r ,70, 
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At ax = 1 the integral in eq. (|79|l is non-singular for £ < 0. In the asymptotic 
regime t — + oo one obtains from (|79|l 

c(i)« C (to)e^ ( '- to) r, (80) 

where the constant value T is equal to 



r^(i- a r s -^Sr /°^^[i-H^« +1 ) . (M, 



For x — > the constant V remains finite for 27 4 ~ 3 > —1, i.e., 27 + 1 > 0. 

The expression 1(801) shows that the expansion rate in Ox 3 direction tends 
asymptotically to the expansion rate in the orthogonal direction, i.e., the uni- 
verse becomes isotropic. The isotropization rate is characterized by the function 
K{t), 

K(t)=logY, *(t) = £W_£W £ (t _>«,)_>(). (82) 
c(t) a(t) 

Since llqi + 4^2 + 93 7^ was assumed in Eq. 1)74(1 . the minimal limit <Zi = 92 = 
q<3 = cannot be taken here, since it implies H a — > 00. To check whether or 
not a corresponding isotropization takes place in the minimally coupled theory 
as well one has to solve the system (|5*5)l - (|5T|l with qi — q 2 = 93 = 0. It is 
straightforward to realize that a solution = A of the minimally coupled 
system requires P(3) = \B 2 . Consequently, a dust equation of state = 
means the absence of the magnetic field. While the non-minimal theory in 
our example admits an isotropization in the presence of a magnetic field, the 
minimally coupled theory does not. 

For a = 1 the initial value problem for c(t) with the initial value c(to) 
degenerates. This requires a special treatment and we do not consider this case 
here. 

6.3.2 "Longitudinal quasi- vacuum": W + P(3) = 

When P( tr )(i) = (7 — l)W(t) and W + P( 3 ) = 0, the conservation law yields 

a(t) 

and a(t) can be found in quadratures from the equation I|tj7|) . The solution 
a(t) = a(to)e Ha ( t ~ t °' > holds in this case as well if the magnetic field term with 
B 2 (t) oc a~ 4 (t) is compensated by the longitudinal pressure P(3)(i) = —W(t) cx 
a~ 2l {t). This requires a transversal stiff matter equation of state, i.e., 7 = 2. 
As in the previous subsection (see the discussion following eq. J73J) we obtain 
two expressions for H a , 

Ha = h ^ H « = 2AP(nq 1+ 4q 2 + q 3 ) > (84) 



W(t) = W(t )(^-) , (S3) 
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which again imply a relation between A and the parameters of the non-minimal 
interaction. From equation (|65|l we find 

c(t) - c(i )e H ° (t ~ to) [l - ae- iH ° <*-*<•)] C [1 - a]" C , (85) 

where 

_ 2^fa)q4 + M 2 (l-2g 1 i? a 2 ) 

C= 8 91 M^2 ' (86) 

Again, the solution for c(t) has the same asymptotical behavior as a(t), i.e., the 
universe becomes isotropic. This solution is non-singular if qi is negative. c(t) 
increases monotonically for 4£ < 1. If 4£ > 1 the function c(t) increases after it 
passed a minimum value. If q\ is positive, there exists a time t* with c(t*) = 0, 
where t* is given by 

t* = t + J-lo g a, a=^. (87) 

The model is then applicable for t > t* . 

Also in this case the isotropization process is a property of the non-minimal 
theory only. The corresponding solution H a = const of the minimal theory is 
necessarily isotropic for W + P($\ = 0. 

6.4 Isotropic universe model with "hidden" magnetic field 

Let us consider now the conditions, under which an isotropic model with a(t) = 
bit) = c(t) is compatible with the existence of a magnetic field. The Einstein 
equations reduce to 

• \ 2 

a 




- =A + KW + KB z it)i- + i9 qi + 4q 2 + q 3 )[-\ }, (88) 



A /-•/',, /••//":/ : J ^ • :>'/! I W.>. I ) I. ''/I % <J.;) 



(89) 



and 



2-+ - =A-KP ( 3 ) +KB J (t)\- + 2 9l --[-\ (13 9 i + 4g 2 + g 3 ) } . (90) 



Furthermore, the trace equation becomes 

2" 



G 



a a 
- + [- 

a \ a 



= 4A+KiW-2P {tl . ) -P (3) )+2 K B 2 it)i9q 1 +4q 2 +q 3 ) 



a la 
a \ a 



(91) 
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and the conservation law simplifies to 

W+^(3W + 2P (tl) +P {3) ) =0. (92) 
Equations (|89|l and l|90|> are compatible when 

In the minimally coupled case compatibility requires Ptgs = P(tr) + -E" 2 - This 
excludes an isotropic matter configuration P/ 3 ) = P(tr) together with a non- 
vanishing magnetic field. For the present non-minimal coupling the situation is 
different. Here, the compatibility condition 193(1 takes the form 

H(t)Q 1+ H 2 (t)Q 2 = l, H(t)=^-, (94) 

a{t) 

where 

Qi = 6(b + 4g 2 + g 3 , Q 2 = 12?i + 4g 2 + q 3 . (95) 

(In the present isotropic case we have iJ a = - = = I = H c = - = H). The 
solutions of this compatibility condition can be classified as follows. 
Case Q\ = 

The Hubble parameter is constant and satisfies the condition 6H 2 qi = 1. 
Case Qi^O 

For positive Q 2 the function H(t) satisfies the relation 



l-VQiH{t) 1-y/QiHito) _Vo2('-'o> 

-e Q i 



1 + VQ^H(t) 1 + VQ^H(t ) 

For the subcase i?(io) = —7=, there exists a special constant solution 



(96) 



H 2 (t)=H 2 (t Q ) = — \ r. (97) 



For negative Q 2 the Hubble function is 



H(t)-±= tan ( arctan (to) + *° } 1 ■ (98) 

For Q 2 = the Hubble parameter is linear in time 

H(t) = H(t ) - . (99) 

Furthermore, for non-minimal coupling the richer structure of the field equations 
admits a configuration for which a non-zero magnetic field does not appear in 
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Einstein's equations. This happens if all the multipliers of the terms nB 2 (t) in 
l|88|) - (|90[) vanish simultaneously. The corresponding conditions are 



H(t) = - = const = H , q x = —^, 4q 2 + q 3 = -—. (100) 

Qj £iS1q 



It is interesting to realize that this solution coincides with l|97|) . The third 
relation admits the particular case 

2 3 

92 = -772 > 93 = 772. ( 101 ) 

for which 6gi + 3q 2 + 53 = and the trace of the susceptibility tensor vanishes 
(see, eq. (JSJ). As a consequence, there exists a stationary cosmological solution 
with 

^(3) = ^(tr) = P -, W(t) = const = W , P = const = - W , 

3Hq = kWq + A , H(t)=H a , a{t) = a(t )e Ha{t ' to) . (102) 

The non- vanishing magnetic field is hidden as far as the space-time evolution is 
concerned. The additional coupling terms give rise to a non-minimal screening 
of the magnetic field. 

Also the solution l|102|) with l|100fl and l|101|) is characterized by a Hubble 
expansion that is directly determined by the non-minimal coupling strength. 
While this solution does not constitute a real inflationary model since there is 
no exit from the de Sitter phase, we hope that it may provide the starting point 
for a more general approach in which the parameters q\ , qi and 53 are no longer 
constants but dynamical degrees of freedom, e.g., a multiplet of scalar fields. In 
such a context more "realistic" solutions might well occur. The circumstance 
that the impact of the non-minimal coupling weakens in the long time limit 
becomes obvious if we try to solve, say eq. I)9()|l . with a power law ansatz a cx t v . 
All the terms in the braces on the right hand side, except for the first one which 
is independent of the non-minimal coupling, decay as t~ 2 . Consequently, these 
terms play a role at early times but they become irrelevant in the long time 
limit. 



7 Discussion 

Inspired by a well motivated non-minimal coupling between gravity and elec- 
tromagnetism we have explicitly demonstrated that the richer structure of the 
corresponding theory gives rise to novel features of the cosmological dynam- 
ics. We have obtained a number of simple exact solutions for Bianchi I models 
with magnetic field. For axially symmetric configurations we found inflationary 
type solutions with magnetic field which describe an isotropization process as 
a result of the non-minimal coupling, i.e., without a counterpart in the mini- 
mally coupled theory. Furthermore, some solutions of the non-minimal theory 
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establish a direct relation between a cosmological constant and the coupling 
parameters of the non-minimal interaction. Finally, we have shown that there 
exists an isotropic de Sitter solution for which the magnetic field is screened by 
the non-minimal coupling. 
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